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We investigate the effect of averaging inhomogeneities on expansion and large-scale structure 
growth observables using the exact and covariant framework of macroscopic gravity (MG). It is 
well-known that applying the Einstein’s equations and spatial averaging do not commute and lead 
to the averaging problem and backreaction terms. For the MG formalism applied to the Friedman- 
Lemaitre-Robertson-Walker (FLRW) metric, the extra term can be encapsulated as an averaging 
density parameter denoted Oyi. An exact isotropic cosmological solution of MG for the flat FLRW 
metric is already known in the literature, we derive here an anisotropic exact solution. Using 
the isotropic solution, we compare the expansion history to current available data of distances to 
supernovae, Baryon Acoustic Oscillations, CMB last scattering surface data, and Hubble constant 
measurements, and find —0.05 < < 0.07 (at the 95% confidence level). For the flat metric 

case this reduces to —0.03 < Da < 0.05. The positive part of the intervals can be rejected if a 
mathematical (and physical) prior is taken into account. We also find that the inclusion of this term 
in the fits can shift the values of the usual cosmological parameters by a few to several percents. 

Next, we derive an equation for the growth rate of large scale structure in MG that includes a term 
due to the averaging and assess its effect on the evolution of the growth compared to that of the 
AGDM concordance model. We find that an term of an amplitude range of [-0.04,-0.02] lead to a 
relative deviation of the growth from that of the AGDM of up to 2-4% at late times. Thus, the shift 
in the growth could be of comparable amplitude to that caused by similar changes in cosmological 
parameters like the dark energy density parameter or its equation of state. The effect could also be 
comparable in amplitude to some systematic effects considered for future surveys. This indicates 
that the averaging term and its possible effect need to be tightly constrained in future precision 
cosmological studies. 

PACS numbers: 98.80.Es,98.80.-k,95.30.Sf 


I. INTRODUCTION 

The rapid improvement in the quantity and quality of incoming and future observational data has encouraged 
the field of astrophysics to aim for a precise and accurate cosmology. In such a context, the understanding and 
control of systematic effects in the data and modeling have become essential to such an endeavor. Indeed, in order 
to constrain the information extracted from the data to a percent precision level, it is necessary to consider nuisance 
effects and contaminants that affect the data at this level. These effects include not only systematics related to data 
measurements but also corrections in the theoretical modeling including nonlinear and relativistic effects. 

One effect worth exploring is the averaging problem in relativity and cosmology that may effect the precision and 
accuracy of the cosmological constraints derived from the data [l|-[2^. The problem originates from the fact that 
spatial averaging (or smoothing out inhomogeneities) in the universe is an operation that does not commute with 
applying the Einstein field equations. In other words, the field equations derived form the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric that describes the universe at very large scales will be different from the equations 
derived at smaller scales and then averaged over large distances and volumes in the universe. This is due to the 
non-linear nature of the field equations of General Relativity. The presence of this non-commutation is well-known 
in the literature and usually gives additional terms in the Friedmann equations known as back-reaction terms 
[13. This back-reaction has been in general agreed upon by now to be too small to affect dramatically the overall 
dynamics of the universe, however it remains an open question whether this back-reaction can affect the constraints 
on cosmological parameters at the percent-precision level and thus if it should be considered at the same footing as 
other systematics in cosmological analyses. In this paper, we address some aspects of this question. 

The averaging procedures and the resulting backreaction terms provide mathematical formalisms on how smaller 
scale inhomogeneities in the universe can affect the dynamics at large scales. Such averaging schemes provide effective 
dynamical equations that explicitly relate the “macroscopic” observables to the underlying “microscopic” structure. 
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The construction of averaging schemes for General Relativity has a long history and there has been many proposed 
metho ds While some have proposed averaging schemes that use only scalar evolution equations, e.g. ii 

MM I14I Il 8 j|. a great deal of effort was put into developing a fully covariant averaging procedure that can be applied 
to tensor equations flan. 

Perhaps one of the most promising covariant schemes was proposed by a a and is referred to as Macroscopic 
Gravity in the literature. This was inspired by some earlier work on the subject matter [3 [ 13 . The formalism 
derives the macroscopic gravitational field equations based on the usual Einstein’s equation plus a tensor term due to 
the averaging process. These new field equations can be solved for a specific macroscopic geometry (for example the 
FLRW metric) without explicit reference to the microscopic geometries. The solutions to the macroscopic equations 
will give the dynamical equations with additional terms due to the averaging that are related to the microscopic 
geometries. This became an attractive framework for cosmology and the Macroscopic Gravity field equations have 
been solved for the flat FLRW metric [3483 some authors have explored the cosmological implications of the 
formalism [ 4 l|, . 

In this paper we analyze the effect of terms due to the macroscopic gravity averaging scheme on constraints from 
the expansion history and the growth rate of large scale structure in the universe. 

The structure of the paper is as follows. In section im we give an overview of the Macroscopic Gravity formalism. 
Then, in section Hill we detail the approach to systematically obtain exact solutions in macroscopic gravity. We briefly 
rederive the flat FLRW isotropic solution and introduce a new anisotropic exact solution. In section El we compare 
the macroscopic gravity observables of the isotropic solution to the expansion history data. In section |V] we derive 
the growth equation for the flat FLRW model. Finally we conclude in section IVTl Units are chosen throughout the 
paper such that c = 1. 


II. THE AVERAGING AND MACROSCOPIC GRAVITY FORMALISM 

In this section we will introduce Zalaletdinov’s Macroscopic Gravity formalism detailed in Refs The Macro¬ 

scopic Gravity formalism consists of a covariant averaging procedure, a method of assigning derivatives to the averaged 
geometric objects, and the application of the averaging procedure to the Einstein Field Equation (EFE) and the Car- 
tan’s Structure Equations in order to derive an effective EFE giving the coarse grained macroscopic dynamics. 

The average of an arbitrary tensor field at a point x over some averaging region Ej, surrounding the supporting 
point x is defined as 


Sx j Tax 

Where, ^/—g{x')d^x' is the 4-volume of the averaging region for support point x, and A°'^,{x, x'), A^ p{x', x) 

are arbitrary bivectors (i.e. two point vectors whose primed index transforms like a vector at x' and unprimed index 
transforms as a vector at x ) which satisfy the conditions \iTa.x'^x A°‘iji{x',x) = and A°‘pi{x,x')A^^„{x',x'') = 
A'^^,,{x,x"). The first condition ensures that the average tensor at a supporting point x (P“^(x)) becomes the value 
of the original tensor field at x when the averaging region goes to zero (limEj,->.o = P°‘p{x)) and the second 

ensures that is the inverse operator of A°^i^{x',x). The most natural and well known bivector is what is 

known as the bivector of geodesic parallel displacement (g“^(x',x)) [13 which satisfies both the required conditions 
and hence can be used as the “averaging bivector”. In fact, when the averaging bivector is the bivector of geodesic 
parallel displacement, the above definition of tensor averaging is similar to the one used by Isaacson in his well known 
paper on gravitational radiation [ 3 , [13 with the exception that in the latter, the integration is over the background 
while in the above definition, the integration is over the actual microscopic geometry, i.e. in the above equation the 
measure y/—g{x') is that of the microscopic metric rather than that of the macroscopic background metric. However, 
when defining the differentiation of average tensors, using the bivector of geodesic parallel displacement introduces 
some complications which we will discuss below, so it won’t be the one used in the Macroscopic Gravity formalism. 

In order to define derivatives of average tensors, the averaging region at each supporting point is defined according 
to the following prescription which Zalaletdinov calls “averaging region coordination”. All the points x°‘ £ Ej, in a 
chosen averaging region (E^,) for a supporting point a;“, are Lie dragged to a nearby supporting point ?/“ = a;“ -I- AAf“ 
along the integral curves of an arbitrary vector field parameterized by A, in order to define the averaging region at 
that point (Ey), using a bivector referred to as the coordination bivector which satisfies the two conditions satisfied by 
the previous averaging bivector. For simplicity, the averaging bivector can be taken as identical to the coordination 
bivector. Now, the averaging region of supporting point y (Ey) read 

Sy = {2/“' |2/“' = + AX^PA'^'pix', x); x'“ £ E,} 
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This procedure is used to construct averaging regions for all the supporting points in the manifold. The well defined 
coordination between the averaging regions allow us to write the measure in a region in terms of the measure in a 
nearby region as 

(l + + 0(AA2) 

Where the semicolon stands for covariant derivation with respect to the microscopic connection. 

Using the above expression, the Lie derivative of the volume can be written as 

=r(Al%)UE. (2) 

With the angle bracket denoting integration over the averaging region divided by Vs, for example =< >. 

The Lie derivative of an arbitrary average vector P“ can be expressed as 

( 3 ) 

where we have defined Pq, = Pa + aP,P' ^ the coma stands for partial derivative, and the over tilde represents the 
bilocal extension of a geometric object, for example P“ = A°‘a'P°‘ ■ 

In order for the average tensors to be single valued functions of the supporting point, the partial derivatives must 
commute. Since £^P°‘ = (P“) ^ — P^^“ ^ the commutator of the partial derivatives is given by 

= (^“:[/ 37 ]) + {P°'-^\'y.Phs) “ [ 7 :/ 3 ]; 5 ) 

Requiring the partial derivatives to commute (P“ = 0) implies that the averaging bivector satisfies the condition 

Furthermore, from Equation ([2]) the condition for the Lie-dragging of an averaging region to be volume preserving 
(i.e. ilfUsx = 0) reads 


A^p.a' = 0 ( 5 ) 

Now, using the conditions (ll]),® and equation ([3]), the partial derivatives of an average vector can be written as 

p^,p = (p'“:^) (6) 

The bivector of parallel propagation does not in general satisfy conditions (|4]) or so it is not possible to set 
up the averaging region coordination using that as the averaging bivector. However, it has been shown [s^ that 
for an arbitrary metric there always exists a bivector satisfying equations ([4]) and m, and that satisfying these two 
conditions is equivalent to the bivector being the product of two vector bases A°‘p,{x,x') = e“-^ (x') with the 

structure functions C\j being constant (where \p(i),e{j)\ = C^ije(k))- Choosing different (i.e. different averaging 
bivectors) will give different averaged fields for a given microscopic tensor field (see Eq. ([I} ). 

In order to obtain an effective EFE, it is necessary to know what averaged geometric object gives the effective 
dynamics. In the Zalaletdinov formalism, the “bilocal extension of the connection coefficients” defined as 



which transforms like a connection at x, like a scalar at x', and reduces to the microscopic connection in the 

limit x' goes to x, is what should be averaged in order to get the effective macroscopic connection coefficient. There 
will be a macroscopic curvature tensor and a macroscopic metric {Gaff) corresponding to the macroscopic 

connection (< >)• Additionally, there exists a connection corresponding to the averaged microscopic 

Riemann tensor The difference between the two connection coefficients is defined as the Affine deformation 

tensor (A“^.^ = < > -7r“^^). 
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By averaging out the Cartan structure equations, the metric compatibility equation , their integrability conditions 
and the microscopic EFE, the Macroscopic Gravity field equations can be constructed and shown to be of the form 




( 8 ) 


Where denotes the macroscopic Ricci tensor, 5 “*^ the average of the inverse microscopic metric, the averaged 
stress energy tensor, and T^arav)'^ ^ gravitational stress energy tensor. 




Here the correlation 2-form is defined as 

7 ^*: ^ _ / -r-a -r-fi 

^ UG X*' £cr] 


-ir 


p[i 




( 9 ) 


( 10 ) 


and it’s traces are as follows. 

Q ol _ o 7€ ol 7e _ o 7€ 

0p(i /3p €7’ ^ p.u'y 

where Q°^pp^ is known as the polarization tensor. 


and Qpi, — Q'^ — Z^ 


The correlation 2-form has the following symmetries 


Z“„, ^ . = 0 

P(7 EP) 


ryOL pL _ rypL OL 

^ (3'y 2/(7 ■^1/7 


ryO /X 

P[l 


= 0 


and it satisfies the equi affine constraint 


^ = 0 
017 2/(7 ^ 


( 11 ) 

( 12 ) 

(13) 

(14) 


The differential properties for the correlation 2-form are set by a correlation 3-form and a correlation 4-form. It is 
possible to set the correlation 3 and 4-forms to zero and hence greatly simplify the formalism by setting the covariant 
derivative of the correlation 2 -form with respect to the macroscopic connection to zero. 


7“ ^ — n 

/3[7 ii<T||A] - ^ 


(15) 


where || represents covariant derivative with respect to the macroscopic connection. This equation also ensures that 
the averaged stress energy tensor is conserved. 

The above equation has the integrability condition. 


“ 7 T\jfe _ n 

P[^l Sktt] - ^ 


7 ^ 7 j\/fa ya 7 , re . ya e 71^7 y 

^ /3[a. Iktt] ~ ^ ^ /3[/i ^ 

Furthermore, setting the correlation 3 and 4 -forms to zero require the quadratic constraint 


(16) 


yo y yO I yO pL y^j Oi _|_ yo 0 7/X tr . yO pL ykj d 

&[l ^ 4e ^ /3[7 kit 5e] P[l va-^ Se KTt] ^[7 2 /(t] 

I yo 6 ryP^ 5 I yo 5 yQ M _ A 

^[1 i^a k-k] ^ /3[7 5e im] ~ ^ 

The average of the Cartan equations implies the Affine deformation tensor needs to satisfy the constraint 


(17) 


4a _ 4 “ — _ln“ 

^ [pa\\p\ ^ e[p^ la] - Ppa 

Additionaly, the average of the integrability condition of the Cartan equations gives 

/3[pl^“e(TA] + /3[pQ“e(TA] “ /3 (tA] “ PaX] = ^ 


(18) 


(19) 


For a given macroscopic metric, equations (un-dill) can be solved to derive the correlation 2 -form and the corre¬ 
sponding additional terms in the field equations. 
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III. EXACT COSMOLOGICAL SOLUTIONS TO MACROSCOPIC GRAVITY EQUATIONS 

In order to solve the Macroscopic Gravity equations, it is essential to assume that the inverse of the averaged 
microscopic metric is equal to the macroscopic metric. This will restrict the class of solutions but it’s not possible to 
avoid this since the theory does not provide a method of deriving this quantity other than explicitly performing the 
averaging. 


A. Algorithmic approach to solving the Macroscopic Gravity equations 

A systematic approach to obtaining an exact solution are as follows. 

• Define the metric for the macroscopic geometry and calculate the Riemannian curvature tensor 

• Define the correlation 2-form in terms of 720 arbitrary functions of the coordinates with the symmetries given 
by equation dm 

• Apply the algebraic cyclic identity Eq. (HI 

• Apply the algebraic equi affine constraint Eq. da 

• Solve the integrability condition Eq. da 

• Solve the differential constraint Eq. da 

• Solve the quadratic algebraic constraint for the correlation 2-form Eq. da 

Solving these equations will in general give all the independent components of the correlation 2 -form (although sym¬ 
metries in the macroscopic geometry can place further constraints reducing the number of independent components). 

• Now the affine deformation tensor can be solved for using equations pSI) and (fa 

• The gravitational stress energy tensor can now be calculate using Eq. (ED 

• Finally any constraints on the gravitational stress energy tensor due to symmetries in the macroscopic geometry 
need to be applied 

All remaining independent functions in the correlation 2-form and the affine deformation tensor will correspond to 
different microscopic geometries giving the same macroscopic geometry, and are free parameters of the model. Now 
the macroscopic EFE (HD can be derived for a given averaged stress energy tensor. 


B. Previously derived spatially homogeneous and isotropic solutions 


The model of cosmological interest is the one with the macroscopic geometry described by the FLRW metric 
in agreement with observations. The macroscopic gravity solutions for the flat FLRW metric has been studied in 
the literature (ssl - fs^ with the second reference giving a systematic analysis of the solutions for the case when the 
correlation 2-form [2 °"and the affine deformation tensor are invariant under the six parameter group 

of Killing vectors (corresponding to the three translational and three rotational symmetries of the metric), and the 
electric part of the correlation tensor is zero. The solution to the correlation 2-form was found to be completely 
specihed by three arbitrary constants A, /12 and bi while the affine deformation tensor was specihed by only A. The 
gravitational stress energy tensor reads: 




/ A! 

a 2 
0 
0 
0 


0 

1 

3 

0 

0 


0 

0 

lA^ 

3 

0 


0 

0 

0 

1 w 

3 


( 20 ) 


We have re-derived this solution in this work and our results are in agreement with the findings of [s^, 1^ . In 
summary, for the macroscopic line element 


ds^ = —dt^ + {dx^ + dy'^ + dz^) 


(21) 
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with an averaged stress energy tensor of the form of a perfect fluid = diag{—p,p,p,p) where p is the energy density 
and p is the anisotropic pressure. 

The macroscopic EFE [8] read 


8ttG 1 a 

3 ^ 3 3 


( 22 ) 


2d d^ 1 A^ 

—+ ^ = -87rGp--4-+A (23) 

a 6 

where over-dots denote partial differentiation with respect to the time coordinate t. 

Hence, the Macroscopic Gravity correlations appear like an extra positive spatial curvature term in the Friedmann’s 
equations. 

The constant A emerges from applying the formalism to the macroscopic flat FLRW metric. It has no explicit scale 
dependence, however, it does implicitly depend on scale in the sense that the derived expression holds only when 
the averaging is performed at a scale large enough to reduce the macroscopic geometry to completely homogeneous 
and isotropic. At smaller scales, the effects due to averaging will not be captured by this simple expression and will 
presumably explicitly depend on the scale. 


C. Spatially homogeneous and anisotropic solutions to Macroscopic Gravity equations 

The exact solutions for Macrosco pic G ravity are known only for the flat homogeneous and isotropic an d the static 
spherically symmetric cases [ssl-l^ Il 22 j| . The non static spherically symmetric solution has been found [l23j | using 
“volume preserving coordinates” and an approximation rather than by solving equations (lll|l - ()19|) directly. In this 
section we will consider the solution for a macroscopically homogeneous, anisotropic and spatially flat metric (i.e a 
macroscopically Bianchi type I metric) of the form 

dS^ = —dt^ + a{t)^dx^ + b{t)^dy^ + c{t)^dz^. (24) 

We note that we will derive this exact solution here just as a further example for macroscopic gravity but we will 
use for the observables and the remaining of the paper the isotropic solution from the previous sub-section. 

We will not assume that correlation 2-form 2°^^^ is invariant under the three parameter group of Killing vectors 
(G 3 ). However, the gravitational stress energy tensor (Eq. ([S])) will be required to be diagonal and invariant under 
the action of G 3 since the average stress energy tensor being considered is invariant under the action of G 3 . In the 
language of [S^l it’s a “Type H” solution. 

The assumptions used to obtain the solutions are as follows, 

• The average of the inverse microscopic metric is equal to inverse macroscopic metric 

= G“^ 

• The averaged microscopic stress energy tensor takes the form 

= diag[-p(t),pi(t),p2(t),P3(t)] 

• The electric part of the correlation 2-form is zero. 

= 0 where = [1, 0 , 0 , 0 ] is the time like vector orthogonal to the hyper-surface of homogeneity 

• The affine deformation tensor will be assumed to be invariant under the action of the group of Killing vectors. 

= 0 where kq) = di and i=x,y,z 

All the following calculations were performed using the publicly available tensor algebra package GRTensor and the 
commercial computer algebra package Maple. 

From equation m the correlation 2-form ostensibly has 720 independent component. We started by defining the 
correlation 2-form with these symmetries in terms of 720 functions of the coordinates. The cyclic identity (1131) gives 
250 independent constraints while the Equi-afiine relation (1141) gives 76 additional independent linear constraints and 
the assumption that the electric part of the correlation tensor is zero gives a further 275, reducing the number of 
independent components to 121. The calculations up-to this point will be true for any metric since the metric and 
it’s derivatives played no role in the equations. So for any macroscopic geometry, the correlation 2-form will have at 
most 121 independent components. 
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Now applying the integrability condition (1161) gives an additional 52 constraints bringing the total number of inde¬ 
pendent components to 69. Solving the differential constraint (fTSl) forces all the functions to be independent of time 
and we are left with 69 functions of the position coordinates. 

The gravitational stress energy tensor can now be calculated and it will be a diagonal and depend on 6 of the func¬ 
tions. Applying the requirement that the gravitational stress energy tensor is invariant under the action of C ?3 gives 
6 differential constraints. Solving them, the gravitational stress energy tensor reads 




1 

a(t)b{t)c{t) 


/ Aa + Bb + Cc 0 
0 Aa 

0 0 

Vo 0 


0 0 \ 
0 0 I 
Bb 0 
0 CcJ 


(25) 


where A, B, C are constants. The correlation 2-form now has 69 independent components, 66 of them are functions 
of the position coordinates and 3 are constants. 

Now the macroscopic EFE ([5]) reads 


ab at be Aa + Bb + Cc ^ 

— H-^-=-;-h SttGp 

ab ac ac abc 


(26) 


b c be A „ „ 

--I-h — ——7- SirGpi 

b c be be 


(27) 


a e ae „ 

-H-1-=-87rGp2 

a c ae ae 


(28) 


a b be C „ ^ 

--b--b — =-r - 87rGp3 

a 0 ab ab 


(29) 


It would be of interest to know how this compares to the dynamics of a homogeneous anisotropic space time with 
spatial curvature. Since the RW metric with positive curvature is a special case of the Bianchi type IX, the natural 
choice for comparison would be of that type. 

The metric for a Bianchi model can in general be written in the form 


dS"^ = —dt^ + dx’^dx'^ 


(30) 


where w^a are the components of the invariant basis I-forms corresponding to the Bianchi type and ’nt pi-i', is a 
symmetric matrix that is a function of only time. For the Bianchi IX model the invariant basis I-forms are [12^ 

= cos{ip)d0 -I- sin(^) sin(d)d(/) 


= sin{'ip)d0 cos(^) sin(0)d(/) 

= dll) + i cos{0)d(j) 

The simplest Bianchi IX model for comparison would be the one with = diag[a{t),bit),e{t)]. For this model 

the EFE read 


db 

ab 

b 

b^ 


at bt 1 
ae ac 


c bt 2a^b^ -I- 2a^c^ -I- 6^ -I- c'’’ — 3a^ — 2b^c^ 
c~^ be Aa'^b'^c^ 


- + 87rGp 

(31) 

— SttGpi 

(32) 


Aa'^b'^c^ 







— 36 '^ — 2a^(? 
4,0^ 


87rGp2 


(33) 


a ^ c ^ ac 


a c ac 


a 

a 



be 

ab 


2a? + 26^ + a'^ + b'^ — 3c^ — 2a^6^ 

4a^62c2 


SirGps 


(34) 


The exact Macroscopic Gravity solutions for the spatially flat anisotropic metric (I26II29I1 have some similarity with 
the spatially closed anisotropic solution for the microscopic EFE. When a{t) = b{t) = c{t), both sets of equations 
have the same form. Equations (I261I29I) reduce to the flat homogeneous solution described in the previous section 
while equations (13111341) reduce to the spatially closed FLRW solution. However, unlike the simple case of the isotropic 
solution, it is unclear how to relate these new terms to a spatial curvatu re. I t is known that the “mixmaster” models 
described by equations (l?T]) - (IM)) show chaotic behavior at early times [l25l| and it remains to be analyzed whether 
the macroscopically anisotropic models have similar behavior that would wipe out any macroscopic anisotropies. 


IV. EXPANSION HISTORY OBSERVABLES AND CONSTRAINTS ON MACROSCOPIC GRAVITY 

ISOTROPIC SOLUTION 

The observational consequence of the macroscopic FLRW model with its additional “dynamical curvature” J? jd? ^ 
on the luminosity distance measurements and hence the constraints on the cosmological parameters from the distance 
observables has been studied in the literature [4l|. In this paper we will apply the additional constraint > 0 
(hence < 0) coming from the constraint equations for the affine deformation tensor (see [1^). We obtain the 
results for the cases with and without this constraint. 

First, we write the Macroscopic RW metric as, 

dS^ = —di? + a{t)^[dr^ + fk{r)^{d9‘^ + sin^ Odej)'^)] (35) 


where 


fk{r) = < 


sin(r) if fc = 1, 
r if fc = 0 , 


and k is the spatial curvature of the macroscopic metric. 


sinh(r) if fc < 0 . 

If the source is a perfect fluid, the effective EFE read: 

= -ttGp - 


A 

"3 


irl! 

3 


(36) 


- = |7rG(p + 3p) + 

a 6 o 


(37) 


where the Hubble parameter has been defined as, H = d/a. The effective Friedman equation (1361) can be rewritten in 
terms of the current matter parameters as 


H{a) = Hoirika-^ + H^a-^ + Ha + H^a"^)" (38) 

where Hm = ^ttGp/Hq is matter density parameter, Ha = A/3Hq is the cosmological constant density parameter, 
Hfc = —kla?Hn IS the curvature density parameter, H^ = —A^/SH^a^ is the “gravitational energy” parameter due 
to averaging [41|, and Hq is the Hubble parameter evaluated today. 

In this work, we make the assumption that light rays on average follow the null geodesics of the averaged macroscopic 
space time and that the only changes to the luminosity distance are due to the change in the modified Friedmann 
equation. Some rays of light are demagnified and some are magnified but on average photon flux conservation leads 
to no net change in the luminosity distance So overall, this seems to be a reasonable assumption for the average 
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of a large number of photons and is consistent with some of the findings in the literature [H, . However some 

authors have found [6i[73 that inhomogeneities could lead to small changes in the redshift relation and possibly 
large changes to the luminosity distance. It is re-assuring that the changes in the redshift were found to be small 
in those studies while the changes in the luminosity distance are suppressed on average by photon flux conservation. 
Nevertheless, it will be good to address this point further in the formalism used in this paper by applying it to the 
null geodesic equation and the Sachs equations explicitly, similar to what was done for the EFE. We leave this full 
project for future work. 

Now, following the usual derivation, the luminosity distance can be written for the curved FLRW macroscopic 
metric as. 


dL = 


1 


-rfk 


^/\^k\da' 


aHoy/l^k 

and for the flat FLRW macroscopic metric it reads 

1 r 


'=1 -b -b + n^a'Y 


dh = 


da' 


aH, 


0 J a' — l 


'=1 -b flAa''^ + rima') 


/4 


(39) 


(40) 


Due the degeneracy between and in the denominator (i.e. Friedmann equation), for the fits to the data 
we use the dynamical energy term that is the sum of the averaging gravitational energy and geometric curvature 
{ilkd = -b Dfc) same as ref (4ll |. 

We fit the cosmological parameters for the FLRW solution of Macroscopic Gravity (and other model^using the 
available cosmological distance data. The supernova (SNe) observations from the Union 2.1 data set [93, Cosmic 
Microwave Background (CMB) last scattering surface data from WMAP 9 year data release [9l|, Baryonic Acoustic 
Oscillations (BAO) from WiggleZ [93 and the Hubble rate from HST measurements 

The parameter hts for the various models were performed using minimization via a maximum likelihood analysis 
(i.e. we minimize = y^^r -b xYad + y^a^fp) and Monte-Carlo Markov Chain approach using a modified version of 
the publicly available package COSMOMC [gj. 


In order to get the constraints from the supernova data, we define xYsn 

2 {^kobs[Zi) - 

XsN — 2 ^- 

1=1 

where = fh — M = Slogj^g [dh^z)) -b 25 is the extinction corrected distance modulus, Ui is the uncertainty in 
the SNe data point, m is the apparent luminosity and is the Luminosity distance measured in Mpc. When 
performing the fits we effectively marginalize over the absolute luminosity M. 

In order to fit for the CMB surface of last scattering we define three fitting parameters [§3 ■ The shift parameter R 
defined by 


R{zY = +z^Da^zY (42) 

the redshift to the surface of last scattering z* given by 

z* = 1048 (1 -b 0.00124(Bt,/i2)-O''^38) (l -b gi{VLmh^fY (43) 

where (see for example [13 ) 


0.0783(Bb/i2)-o-238 
1-b 39.5(Bf,/l2)0.763 

0.560 

l-b21.1(Bb/i2)i-8i 


and the acoustic scale (Zq) defined as 


I 


a 


(1 -b 2») 


t^Da{z*) 

rs{z*) 


(44) 
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with the proper angular diameter distance, Da{z) = dL{z)/{\ + and the comoving sound horizon 

If TTJ7 da 

rs{z*) = 


V 3 J 0 a^H{a)^yi + (30b/417^)a 


where = 2.469 * for Tcmb = 2.725K. 

The parameters are fitted using 


XCMB=^^iC0V ^{Xi,Xj)Axj 


(45) 


(46) 


with Xi = (R, Ig, z^,), Axi = Xi — x^ and Cov~^{xi,Xj) the inverse covariance matrix for the parameters from the 
reference 

Next, in order to obtain the constraints from the BAO, following we define the effective distance Dy as 


D,iz)=[Dliz)il + z) 


with the redshift at the decoupling epoch given by 

1291(n™h2)0.25i 


Zd = 


H{z) 

(l + 6 i(Ob/i 2 )) 


with 


1 + 0.659(n™/i2)0-828 

hi = 0.313(0^h2)-o-4i9 (1 + 0.607(fl™/i2)O-674) 


(47) 


(48) 


62 = 0.238(fl^/l2)0'233 

The parameter constraints from the BAO are now given by 


Xbao 


= E 


si^Zi^ j D y{^Zi) {x si^Zi^ ! Dyi^Zjfj^Qds 


(49) 


We also add the prior, the inverse of the angular diameter distance at red shift 0.04 equals 6.49405 x lO”^ ± 
0.31512 X 10“9^ (i.e. Hq = 74.2 ±3.6 km/s for the fiducial model) given by the HST measurements [o^, and the prior 
Q.bh? = 0.022 ± 0.002 given by big bang nucleosynthesis. 

We perform the parameter hts for Macroscopic Gravity, ACDM and wCDM with a constant equation of state, by 
varying the physical dark matter density (UuMh^), the physical baryon density (Vldh?), the curvature parameter (Ofc) 
and in the cases of Macroscopic Gravity and wCDM, the dynamic curvature i.e. the sum of the averaging gravitational 
energy and geometric curvature (fifed = ± fife) and the equation of state of dark energy (w) respectively. The 

values for those parameters and the derived parameters flm, ^A^ f^A and Hq are summarized in TableHl and Figure[T] 

For the macroscopic gravity, we find in the case of models restricted by the mathematical and physical prior 
[S^ . that —0.027 < ± 0 (68% confidence level). In the case where we do not impose the prior, we obtain 

—0.024 < flA < 0.036. As we will discuss further in section EH the mathematical prior turned out also to be 
a physical prior consistent with the fact that a larger magnitude of negative backreaction term leads to a larger 
enhancement of the growth of structure as supported by studies using inhomogeneous cosmological models [9^, . 

We are also able to reproduce exactly the results of reference lo w here the table shows that the constrained value 
for Ha is significantly large when the SDSS SNelA compilation |l0Cll | (which uses the MLCS2k2 light curve fitter) is 
used. TableHl and Figure [T] uses only SNe data from the Union 2.1 compilation that uses the SALT II light curve 
fitter. 


V. GROWTH OF LARGE SCALE STRUCTURE IN THE MACROSCOPIC GRAVITY AVERAGED 

UNIVERSE 

A. Derivation of the growth evolution equations in Macroscopic Gravity 

In order to study the growth of structure due to small inhomogeneities in the macroscopically flat Friedmann 
universe, we will perturb the metric, the stress energy tensor, the correlation 2-form and the affine deformation tensor 
about the exact macroscopic solution. The new quantities will be given by 

9°‘f3 =9^°^°‘B+Sg' 


(50) 
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Parameters 

ACDM 

wCDM 

MG with prior 

MG w/o prior 

Flat MG with prior 

Flat MG w/0 prior 

Hk 

0.014«“« 

-0.0i3l°:°)7 

U.UUZO_Q QQ4y 

U.UU( O_0.0057 

0 

0 

Qa 

0.6881°°" 

0.690l°:°l° 

0.7241°;°?° 

U.DyZ_o.036 

n 7QQ+0-016 

U. / OO_0 

n 71 1 +0-023 

U. / -l--l-_0.023 


0.298l°°f5 

0.297l°;°g 

0.295l°;°ll 

u.zi;o_Q 092 

0.279l°;°l° 

0.280l°;°l° 

Ho 

71.51=1 

72.8ll,° 

69.9l?;^ 

69-9ll:o 

69.5111 

72.8ll,° 


N/A 

N/A 

—0 021 fi+0.0216 

U.UZiO_o.o054 

0.0058l°:°i° 

U.Uizo_o 0093 

0.009l°;°l° 

w 

-1 

— 1 

-1 

-1 

-1 

-1 


TABLE 1: Marginalized parameter constraints (68 % confidence) from the cosmological distance observations (supernova data 
from the Union2.1 compilation, the HST data, the last scattering surface data from WMAP9, and the WiggleZ BAO data). 
The prior here is flA < 0.. Results are given for the Macroscopic Gravity using the curved and fiat macroscopic FLRW metric, 
with and without the prior. 


rpct _ 


5T° 


(51) 


ya _ 7 ( 0 )“ 

^ UG ^ /37 yo 


0-i va 


(52) 




(53) 


where the superscript zero denotes the value of the exact solution and the prefix <5 denotes the perturbations about 
the exact value. The perturbations are taken to be of order e . By expanding the macroscopic gravity equations 
(|8 |) - (I19I) in terms of the order parameter up to linear order, it’s possible to obtain the equations governing the first 
order terms. The perturbations we consider will be on scales smaller than the averaging domain, so they will be the 
fluctuations that would be smoothed out from the averaging. 

We reproduced the perturbation results of [s^, but present here the derivation steps, and we derive the growth 
rate equation further below. 

The zeroth order terms will satis fy th e original exact equations while the first order terms will satisfy a linearized 
version of the above equations (see |10l| ). The equations (ITT]) - (IT9l) at first order will be given by 


^ =0 
/3(7 ycr) 


(54) 


/i _ TyfJ, a 
/37 1/(7 — 1/7 Pa 


P[l 


= 0 


=0 
(3(7 va ^ 


(55) 

(56) 

(57) 


,,, =0 
^[7 yG\X] 


/ 3 [fi qy ckttJ e[^ ph/TrJ p[^ dy ±k, 7 t] p[fi ey dK 


Zio) 


e 7 


6M° 




- 


a 7 
e[p Sy 






(58) 

(59) 


(6/ S 


S 7(0)“ ^ 1^ 7(0)^ “ ^ 7(0)^ ^ _i_A 7 “ ^ 7(0)^ 

/3[7 Se 1/(7] + 0^ /3[7 ua'^ ktt £e] + 0^ /3[7 1 /( 7 ^ Se k-w] + p[-y Se k-w ua\ 


ya 

-iM S 


K7Z 8e 




-,(5 /7 
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0.1 26 
0.1 24 
0.122 
0.1 2 
0.1 1 8 
0.1 1 6 








FIG. 1: Two dimensional marginalized joint contour plots ((68% and 95% confidence levels)) for the FLRW solution of the 
Macroscopic Gravity using the supernova data from the Union2.1 compilation, the HST data, the last scattering surface data 
from WMAP9, and the WiggleZ BAG data. These are the results for the Macroscopic Gravity model without any prior on 
and where spatial curvature is also allowed. 


(61) 

0[pQ\.x] - 

where | represents covariant derivative with respect to the unperturbed macroscopic connection. 

If Z^°'> ~ 0(e) (i.e. J\?,hi,b 2 ~ 0(e)) the equations governing the first order correlation 2-form and 

the affine deformation tensor will be identical to the ones satisfied by the zeroth order quantities CID-dlll). Therefore 
if we assume that electric part of the first order correlation tensor is zero = 0 where v'^ = ^(1 — (j),Su'^) ) 

the first order gravitational stress energy tensor will have the form of a positive spatial curvature term. 

The assumption that ~ 0(e) is consistent with the small mean value for n _4 obtained in section HVl 

Since only scalar perturbations are relevant for the growth of inhomogeneities we will restrict the metric perturba¬ 
tions to just the scalar part. Now, without any loss of generality the metric can be written in the conformal Newtonian 
gauge as 

dS'^ = a(? 7 )^(—(1 -I- 2(j))dri'^ + (1 ~ 2ij;){dx'^ + dy^ + dz^)) (63) 

The source will be considered as a perturbed perfect fluid and the first order stress energy tensor will read 

5f\ = -Sp 
= ^ip + p)Su, 
dTG = SpS) 


(64) 

(65) 

( 66 ) 
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where 5p is the energy density perturbation, dp is the pressure perturbation and Sui is the comoving peculiar 
velocity 

Now the modified Einstein field equation ([8|) at first order read 


- 3n{ncP + (P') = AirGa^idp + dpA) 

(67) 

yi{U(p + (p') = -AnGa^ [p + P - j ^ idu 

(68) 

(p" + m(p' + { 2 n' + n'^)(p = AnGa^ (dp - 

(69) 


where is the spatial covariant derivative, Spy^ is the energy perturbation to the gravitational stress energy tensor, 
a prime denotes the derivative with respect to ry, is the irrotational part of the comoving peculiar velocity of 

the fluid (which can be written as a divergence of a function Su) and H is defined as Equation (l67ll is the (0,0) 
component, equation (1551) is the (i,0) component and equation (1551) is the (i,j) component where iy^ j . When there is 
no anisotropic stress the i = j component gives ip = (j) and has been used to eliminate ip from the above equations. 

At this point it is convenient to decompose the perturbations into the eigenfunctions of the Laplace equation (see 
[l02| '). Scalar harmonics satisfy + k^Q = 0 while vector harmonics are given by Qi = . 

It becomes apparent that for subhorizon modes (k ?7 »1) (p and (p' are negligible compared to the spatial derivatives 
of (p. Hence equation (1671) rewritten as 


W^(P = ATTGa'^p(^S+^'^ 


The twice contracted Bianchi identity for the modified EFE gives 


rpa , rp(grav)a 

^ /3||“ /3||ct 


= 0 


(70) 


(71) 


It might at first appear that the stress energy components are not independently conserved. However that is not the 
case. It can be shown that the differential constraint on the correlation 2 form (USD implies the gravitational stress 
energy tensor is conserved {T^^rav) hence, the averaged stress energy tensor is conserved independently. 

The first order conservation equations for the averaged stress energy tensor are given by 

dp' + 3'H{dp + dp) — 3(p'{p + p) + a{p + p)V^du = 0 (72) 

4((p + p)a'^S/‘^du)’ + V^dp +{p + p)V^(P = 0 (73) 


Where, the first equation comes from <5T'“p|^ = 0 and the second equation comes from the spatial divergence of 
dT" 1 ^ = 0. Since we are considering the matter dominated era, the radiation can be neglected. Hence, p = 0, and 
(po^) is a constant. Defining the density contrast by dm = dp/p and using equations (1721) and (l73)) . we can obtain an 
obtain an evolution equation for the density contrast of the form: 


-nd„ 


AttGo? p 



= 0 


(74) 


In order to proceed, it is necessary to write the perturbation to the gravitational energy density in terms of the 
matter energy density. In order to do that we argue that even though the matter stress energy tensor and the 
gravitational stress energy tensor are conserved independently, the perturbation to the gravitational energy density 
must be tightly coupled to the perturbation of the matter energy density, since the inhomogeneities in the matter 
cause the gravitational stress energy. Hence, we can assume that the comoving peculiar velocity of the gravitational 
energy density is the same as the matter comoving peculiar velocity. With this assumption, the = 0 

component of the first order conservation equations of the gravitational stress energy tensor is given by 

{dA - 2(py + ^a\7du = 0 (75) 

O 


where 


<5^ = 


SpA 


L SttG J 


(76) 
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FIG. 2: The relative difference in growth factor between the Macroscopic Gravity model MG-ACDM and the AGDM. A non¬ 
zero contribution of the Qa backreaction term due to averaging increases the growth proportionally to the amplitude of this 
term. The effect of is up to 2-4% at late times on the growth factor function compared to the AGDM. 


Using equations (l75)l and (|7^ to eliminate Su and integrating, since Ja = 0 when Sm = 0, we find: 

- '^5A = 0 

Which gives 

c '-’m 

~ ~ UttGo^ 

Substituting the above in equation (TM)) gives the growth equation 


- ^47rGa^p - 5m = 0 

The above equation in terms of the cosmological time t, reads: 

/12 \ 

Where the dot denotes differentiation with respect to t, and H is the Hubble parameter d/a. 

In order to conveniently plot the growth, the growth equation can be written as a function of scale factor 

C + (( 1 » )' + 0 = 0 

where prime now denotes partial differentiation with respect to a rather than rj. 


Sm + 2H5„ 


— [AnGp 


(77) 


(78) 


(79) 


B. Effect of averaging on the growth rate versus precision cosmology requirements 


In order to compare the growth of structure within a given model to observational data, it is most common to use 
the log arithmic g rowth factor since that is what is measured from for example redshift distortions and Lyman-Alpha 


forests [iMIml 


/ = 


d In (5 


d In a 


(80) 
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FIG. 3: LEFT; The relative difference in the growth factor function between the ACDM model with Ha = 0.70 and those where 
Ha takes the values shown on the figure. Increasing Ha increases the late-time growth suppression as expected. RIGHT: The 
relative difference of the growth factor between the ACDM model (i.e. the equation of state ui = — 1) and dark energy models 
where the equation of state takes the values shown on the figure. In both cases, the growth factor function is changed by up 
to several percent. 


It would be of interest to know how the effects of back-reaction on the growth compares to that of the dark energy 
density of equation of state. Particularly whether the effects of averaging can be degenerate with a c hange in these 
parameters. For dynamical dark energy models with a constant equation of state w (see for example (ll2Ml2nl| L the 
growth equation is given by 


df 

d In a 


+ f + 





(81) 


with 




= 


,-3 




§ 2 = 1 ^ (-3F!^a-4 - 211^0-3 - 3(1 + u;)OAa-3(i+-)) 

Substituting <5 = f / and S = ^(/^ ~ / + Jin a ') ™ equation (17^ we obtain the growth equation for the Macroscopic 
Gravity model in terms of the growth factor as 


where 


df 

d In a 


+ f + 





-t“ 



^ = 0 


-^0 


(82) 


^ Sflma — 2fliya ^ — 2fl^a 

We find that a non-zero negative term of 2 — 4% due to averaging has the effect of enhancing the growth by 
2 — 4% at late times relative to when no averaging backreaction is taken into account (see Fig. [T]). These effects are 
of the same order as those resulting from changing the dark energy density parameter or its equation of state (see 

Fig. [SD. 
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VI. CONCLUSION 

In this work we studied the effects of averaging inhomogeneities on the expansion history and the growth rate of 
large scale structure using the non perturbative framework of Macroscopic Gravity. The framework is based an exact 
mathematical formalism developed to provide a covariant averaging procedure. The formalism results in modified 
Friedmann equations with a new term that can be viewed as a back-reaction term and have been previously called as 
the averaging gravitational energy density parameter 17^. 

As examples of exact solutions to macroscopic gravity held equations, we rederive here a previous isotropic and 
homogeneous solution and we obtain a new homogeneous but anisotropic solution. Starting from the macroscopically 
homogeneous, anisotropic, and spatially hat metric of Bianchi type-I, we derive the effective Einstein held equations 
with new terms due to the averaging process. These dynamical equations have the form of an anisotropic generalization 
to the Friedmann equations obtained for the isotropic solution and reduce to them when isotropy of the scale factor 
is restored. Unlike the simple case of the isotropic solution, it is unclear how to relate these new terms to a spatial 
curvature. We use for comparison to observations the isotropic solution. 

We then compare the macroscopic gravity expansion equations to available data sets from distances to supernovae, 
Baryon Acoustic Oscillations, CMB last scattering surface data, and Hubble constant measurements. We note that 
instead of using directly Ha in the analysis, we use a term that takes into account its degeneracy with spatial curvature 
into the Friedmann equation so Ha become a derived parameter. 

We hnd for the isotropic macroscopic FLRW metric solution —0.05 < Ha < 0.07 (at the 95% conhdence level). In 
the hat metric case, the bounds reduce to —0.03 < Ha < 0.05. If we take into account a mathematical and physical 
prior that restricts the sign of the averaging term to be negative, then the positive part of the interval can be rejected 
leading to tighter constraints. It is worth noting that the other cosmological parameters (Ha, Hm, Hk, and Hq) are 
moved by a few to several percent from their ACDM concordance model values wheu the averaging term is included 
in the analysis. 

Next, we explore the effect of the averaging term on the growth rate of large scale structure. We rederive previous 
results from perturbing the Macroscopic Gravity held equations and then derive a growth rate equation that can be 
compared to future observations. We assess the effect of the amplitude of the resulting averaging term on the growth 
rate function and hnd that an Ha term of of amplitude range interval [-0.04,-0.02] lead to an enhancement deviation 
of the growth up to 2-4% at late times. This change in the growth is comparable in amplitude to the changes that 
will be caused by a similar change in the dark energy density parameter or its equation of state. Particularly, the 
effect of increasing the magnitude of the negative averaging term is to enhance the growth rate of large scale structure 
which is physically consistent with o ther results in the literature studying gravitational infall/clustering using exact 
inhomogeneous cosmological models [l2l| . The effect of the averaging on the growth is also comparable in amplitude 
to some systematic affects in ongoing and future surveys. 

We conclude from using the averaging macroscopic gravity formalism to assess the effect of inhomogeneities on the 
expansion history and the growth rate of structure that this effect needs to be tightly constrained and analyzed in 
the future for a precise and accurate cosmology. 


Acknowledgments 

We thank J. Dossett for useful comments on modifications to GAMB and GosmoMG. MI acknowledges that this 
material is based upon work supported in part by NASA under grant NNX09AJ55G and an award from the John 
Templeton Foundation. Part of the calculations for this work have been performed on the Gosmology Gomputer 
Gluster funded by the Hoblitzelle Foundation. 


[1] R. Zalaletdinov, Gen. Rel. Grav. 24 (1992), 1015-103. 

[2] R. Zalaletdinov, Gen. Rel Grav. 25 (1993), 673-695. 

[3] J. Behrend, (2008), PhD Thesis, arXiv:0812.2859 [gr-qc]]. 

[4] J. Boersma, Phys. Rev. D.57, (1998), 798-810. 

[5] J. Brannlund, R. v. d. Hoogen and A. Coley, Int. J. Mod. Phys. D 19, 1915 (2010) arXiv:1003.2014 [gr-qc]]. 

[6] T. Buchert and J. Ehlers, Astron. Astrophys. 320 (1997), 1-7. 

[7] T. Buchert and M. Garfora, Phys. Rev. Lett. 90 (2003), 031101 (1-4). 

[8] T. Buchert, Gen. Rel. Grav. 32 (2000), 105-125. 

[9] T. Buchert, Gen. Rel. Grav. 33 (2001), 1381-1405. 

[10] M. Garfora and M. Piotrkowska, Phys. Rev. D.52 (1995), 4393-4424. 



17 


[11] A. A. Coley, Class. Quant. Grav. 27 , 245017 (2010) arXiv:0908.4281 [gr-qc]]. 

[12] F. Debbasch, Eur. Phys. J. B, 37(2), (2004), 257. 

[13] T. Futamase, Phys. Rev. Lett. 61 (1988), 2175-2178. 

[14] T. Futamase, Phys. Rev. D.53 (1996), 681-689. 

[15] M. Gasperini, G. Marozzi, F. Nugier and G. Veneziano, J. Gosmol. Astropart. Phys. 07 (2011), 008. 

[16] R. Isaacson, Phys. Rev. 166 (1968), 1263-1271. 

[17] R. Isaacson, Phys. Rev. 166 (1968), 1271-1280. 

[18] M. Kasai, Phys. Rev. D.47 (1993), 3214-3221. 

[19] M. Korzynski, Glass. Quantum Grav. 27 (2010), 105015 (21pp). 

[20] T. Noonan, Gen. Rel. Grav 16 11 (1984), 1103-1118. 

[21] T. Noonan, Gen. Rel. Grav 17 6 (1985), 535-544. 

[22] A. Paranjape and T. Singh, Phys. Rev. D.76 (2007), 044006. 

[23] A. Paranjape, Phys. Rev. D.78 (2008), 063522. 

[24] M. Shirokov and I. Fisher, Sov. Astron. J. 6 (1963), 699-705. 

[25] N. Zoltov and R. Stoeger Class. Quantum Grav. 9 (1992), 1023-1031. 

[26] N. Zoltov and R. Stoeger, Astrophys. J. 453 574. 

[27] Baumann, D., Nicolis, A., Senatore, L., & Zaldarriaga, M. 2012, J. Cosmol. Astropart. Phys., 7, 051 arXiv: 1004.2488 
[astro-ph.CO]]. 

[28] L. R. W. Abramo, R. H. Brandenberger and V. F. Mukhanov, Phys. Rev. D 56, 3248 (1997) gr-qc/9704037 . 

[29] S. Andrianomena, C. Clarkson, P. Patel, O. Umeh and J. P. Uzan, JCAP 1406, 023 (2014) arXiv: 1402.4350 [gr-qc]]. 

[30] Bagheri, S., & Schwarz, D. J. 2014, J. Cosmol. Astropart. Phys., 10, 073 

[31] J. Behrend, I. A. Brown and G. Robbers, JCAP 0801 , 013 (2008) arXiv:0710.4964 [astro-ph]]. 

[32] D. Bertacca, R. Maartens and C. Clarckson, (2014), astro-ph.CO 1405.4403 

[33] S. Bildhauer and T. Futamase, Gen. Rel. Grav. 23 (1991), 1251-1264 

[34] C. Boehm and S. Rasanen, JCAP 1309, 003 (2013) arXiv:1305.7139 [astro-ph.CO]]. 

[35] K. Bolejko and R. A. Sussman, Phys. Lett. B 697, 265 (2011) arXiv: 1008.3420 [astro-ph.CO]]. 

[36] 1. A. Brown, A. A. Coley, D. L. Herman and J. Latta, Phys. Rev. D 88, 083523 (2013) arXiv: 1308.5072 [gr-qc]]. 

[37] 1. A. Brown, J. Latta and A. Coley, Phys. Rev. D 87, no. 4, 043518 (2013) arXiv:1211.0802 [gr-qc]]. 

[38] P. Bull and T. Clifton, Phys. Rev. D 85, 103512 (2012) arXiv:1203.4479 [astro-ph.CO]]. 

[39] C. Clarkson and O. Umeh, Class. Quant. Grav. 28, 164010 (2011) arXiv:1105.1886 [astro-ph.CO]]. 

[40] C. Clarkson, G. F. R. Ellis, A. Faltenbacher, R. Maartens, O. Umeh and J. P. Uzan, Mon. Not. Roy. Astron. Soc. 426, 
1121 (2012) arXiv: 1109.2484 [astro-ph.CO]]. 

[41] C. Clarkson, T. Clifton, A. Coley and R. Sung, Phys. Rev. D 85, 043506 (2012) arXiv:1111.2214 [astro-ph.CO]]. 

[42] Clifton, T., Gregoris, D., Rosquist, K., & Tavakol, R. 2013, J. Gosmol. Astropart. Phys., 11, 010 

[43] Di Dio, E., Vonlanthen, M., & Durrer, R. 2012, J. Cosmol. Astropart. Phys., 2, 036 

[44] G. F. R. Ellis, General Relativity and Gravitation, ed. B Bertotti, F de Felice and A Pascolini (Dordrecht: Reidel, 1984), 
215-288 

[45] K. T. Enqvist and T. Mattsson, JGAP 02 (2007) 019. 

[46] E. E. Flanagan, Phys. Rev. D 71, 103521 (2005) hep-th/0503202 . 

[47] G. Geshnizjani, D. J. H. Ghung and N. Afshordi, Phys. Rev. D 72, 023517 (2005) astro-ph/0503553 . 

[48] S. R. Green and R. M. Wald, Phys. Rev. D.83 (2011), 084020 (1-27). 

[49] G. Hellaby, Gen. Rel. Grav. 20 (1988), 1203-1217 

[50] G. M. Hirata and U. Seljak, Phys. Rev. D 72, 083501 (2005) astro-ph/0503582 . 

[51] A. Ishibashi and R. M. Wald, Class. Quantum Grav. 23 (2006), 235-250. 

[52] M. Kasai, Prog.Theor.Phys. 117 (2007) 1067-1075. 

[53] E. W. Kolb, S. Matarrese and A. Riotto, New J. Phys. 8 (2006), 322 (1-25) 

[54] J. Larena, Phys. Rev. D.79 (2009), 084006 (1-6) 

[55] J. Larena, J. M. Alimi, T. Buchert, M. Kunz and P. S. Gorasaniti, Phys. Rev. D 79 (2009), 083011 (1-15). 

[56] M. Lavinto, S. Rsnen and S. J. Szybka, JCAP 1312 , 051 (2013) arXiv:1308.6731 [astro-ph.CO]]. 

[57] N. Li and D. J. Schwarz, Phys. Rev. D 78, 083531 (2008) arXiv:b710.5073 [astro-ph]]. 

[58] V. E. Marra, E. W. Kolb and S. Matarrese, Phys. Rev. D 77 (2008), 023003 (1-13). 

[59] P. Martineau and R. Brandenberger, astro-ph/0510523 

[60] T. Mattsson and M. Mattsson, JCAP 0802 , 004 (2008) arXiv:0708.3673 [astro-ph]]. 

[61] G. Montani, R. Ruffini and R. Zalaletdinov, Glass. Quant. Grav. 20 , 4195 (2003) gr-qc/0307077 . 

[62] Y. Nambu and M. Tanimoto, (2005), gr-qc 0507057vl 

[63] F. Nugier, (2013), arXiv:1309.6542r^l 

[64] A. Notari, (2005), astro-ph/0503715 

[65] U. Obinna, J. Larena and G. Glarckson, JCAP 1103 (2011) 029 

[66] A. Paranjape and T. P. Singh, Phys. Rev. Lett. 101 , 181101 (2008) arXiv:0806.3497 [astro-ph]]. 

[67] P. Premadi, H. Martel, R. Matzner and T. Futamase, Astrophys. J. Suppl. 135, 7 (2001) astro-ph/0101359 . 

[68] S. Rasanen, JCAP 0804, 026 (2008) arXiv:0801.2692 [astro-ph]]. 

[69] S. Rasanen, J. Cosmol. Astropart. Phys. 03 (2010), 018. 

[70] S. Rasanen, 2012, Phys. Rev. D , 85, 083528 

[71] M. Reiris, Class. Quantum Grav. 25 (2008), 085001 (26pp). 


18 


[72] S. Weinberg, Astrophys. J. 208, LI (1976). 

[73] H. Russ, M. H. Soffel, M. Kasai and G. Borner, Phys. Rev. D 56, 2044 (1997) astro-ph/9612218 . 

[74] J. D. Schwarz, arXiv;1003.3026 [astro-ph.CO] 

[75] U. Seljak and L. Hui ASP Conference Series vol 88, p 267 

[76] M. Serin, Class. Quantum Grav. 18 (2001), 5329-5352. 

[77] E. R. Siegel and J. N. Fry, Astrophys. J. 628, LI (2005) astro-ph/0504421 . 

[78] R. A. Sussman, Class. Quant. Grav. 28, 235002 (2011) arXiv:1102.2663 [gr-qc]]. 

[79] H. Tanaka and T. Futamase, Prog. Theor. Phys. 117, 183 (2007) astro-ph/0612151 . 

[80] M. Tanimoto, Prog. Theor. Phys. 102 (1999), 1001. 

[81] C. Wetterich, Phys. Rev. D 67, 043513 (2003) astro-ph/0111166 . 

[82] A. Wiegand and T. Buchert, Phys. Rev. D 82 (2010), 023523 (1-24). 

[83] D. L. Wiltshire, New J. Phys. 9 (2007), 377 (1-66). 

[84] M. Mars and R. M. Zalaletdinov, J. Math. Phys. 38 (1997) 4741 dg-ga/9703002 ; 

[85] A. A. Coley, N. Pelavas and R. M. Zalaletdinov, Phys. Rev. Lett. 95 (2005) 151102 gr-qc/0504115 . 

[86] R. J. van den Hoogen, J. Math. Phys. 50 (2009) 082503 arXiv:0909.0070 [gr-qc]]. 

[87] T. Clifton, A. Coley and R. V. D. Hoogen, JCAP 1210, 044 (2012) arXiv:1209.1085 [astro-ph.CO]]. 

[88] J. L. Synge, North-Holland, Amsterdam, 1960 

[89] Weinberg, S. 1976, Astrophys. J. Lett., 208, LI 

[90] N. Suzuki, D. Rubin, C. Lidman, G. Aldering, R. Amanullah, K. Barbary, L. F. Barrientos and J. Botyanszki et ai, 
Astrophys. J. 746, 85 (2012) arXiv:1105.3470 [astro-ph.CO]]. 

[91] G. Hinshaw et al. [WMAP Collaboration], Astrophys. J. Suppl. 208 (2013) 19 arXiv:1212.5226 [astro-ph.CO]]. 

[92] C. Blake, E. Kazin, F. Beutler, T. Davis, D. Parkinson, S. Brough, M. Colless and C. Contreras et al., Mon. Not. Roy. 
Astron. Soc. 418 (2011) 1707 arXiv:1108.2635 [astro-ph.CO]]. 

[93] A. G. Riess, L. Macri, S. Casertano, M. Sosey, H. Lampeitl, H. C. Ferguson, A. V. Filippenko and S. W. Jha et al., 
Astrophys. J. 699 (2009) 539 arXiv:0905.0695 [astro-ph.CO]]. 

[94] A. Lewis and S. Bridle, Phys. Rev. D 66 (2002) 103511 astro-ph/0205436 . 

[95] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) arXiv:1001.4538 [astro-ph.CO]]. 

[96] W. Hu and N. Sugiyama, Astrophys. J. 471, 542 (1996) astro-ph/9510117 . 

[97] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. J. 633 (2005) 560 astro-ph/0501171 . 

[98] A. Peel, M. Ishak and M. A. Troxel, Phys. Rev. D 86 (2012) 123508 arXiv:1212.2298 [astro-ph.CO]]. 

[99] K. Bolejko, Phys. Rev. D 75 (2007) 043508. 

[100] R. Kessler, A. Becker, D. Cinabro, J. Vanderplas, J. A. Frieman, J. Marriner, T. MDavis and B. Dilday et al., Astrophys. 
J. Suppl. 185, 32 (2009) arXiv:0908.4274 [astro-ph.CO]]. 

[101] R. M. Wald, Chicago, Usa: Univ. Pr. (1984) 491p 

[102] H. Kodama, M. Sasaki, Prof. of. Theor. Phys. Sup. 78 (1984) 

[103] L. Guzzo et al.. Nature, 451, 541 (2008). 

[104] M. Colless et al., Mon. Not. R. Astron. Soc. 328, 1039 (2001). 

[105] C. Blake et al., Mon. Not. R. Astron. Soc. 415, 2876 (2011). 

[106] M. Tegmark et al., Phys. Rev. D 74, 123507 (2006). 

[107] N. P. Ross et al., Mon. Not. R. Astron. Soc. 381, 573 (2007). 

[108] J. da Angela et at, Mon. Not. R. Astron. Soc. 383, 565 (2008). 

[109] M. Viel, M. G. Haehnelt, and V. Springel, Mon. Not. R. Astron. Soc. 354, 684 (2004). 

[110] M. Viel and M. G. Haehnelt, Mon. Not. R. Astron. Soc. 365, 231 (2006). 

[111] P. McDonald et ai, Astrophys. J. 635, 761 (2005). 

[112] J. Dossett, M Ishak, Y Gong and A. Wang, J. Cosmol. Astropart. Phys. 04 (2010) 022. 

[113] E.V. Linder, Phys. Rev. D 72, 043529 (2005). 

[114] D. Polarski and R. Gannouji, Phys. Lett. B 660, 439 (2008). 

[115] Y. Gong, M. Ishak, A. Wang, arXiv:0903.0001vl [astro-ph.CO] (2009). 

[116] M. J. Mortonson, W. Hu and D. Huterer, Phys. Rev. D 79, 023004 (2009). 

[117] E.V. Linder and R.N. Cahn, Astropart. Phys. 28, 481 (2007). 

[118] R. Gannouji and D. Polarski, J. Cosmol. Astropart. Phys. 05 (2008) 018. 

[119] Y. Gong, M. Ishak, and A. Wang, Phys. Rev. D 80, 023002 (2009). 

[120] M. Ishak, J. Dossett, Phys. Rev. D 80, 043004 (2009) 

[121] M. A. Troxel, A. Peel, M. Ishak, J. Cosmol. Astropart. Phys. 12 (2013) 048. 

[122] R. J. Van Den Hoogen, Gen. Rel. Grav. 40, 2213 (2008) arXiv:0710.1823 [gr-qc]]. 

[123] A. A. Coley and N. Pelavas, Phys. Rev. D 74, 087301 (2006) astro-ph/0606535 . A. A. Coley and N. Pelavas, Phys. Rev. 
D 75, 043506 (2007) gr-qc/0607079 . 

[124] H. Stephani, E. Herlt, M. MacCullum, C. Hoenselaers and D. Kramer, “Exact Solutions of Einstein’s Equations,”; 
L. D. Landau and E. M. Lifschits, “The Classical Theory of Fields : Course of Theoretical Physics, Volume 2,”. 

[125] C. W. Misner, Phys. Rev. Lett. 19, 533 (1967). 





